Spontaneous emission enhancement in metal-dielectric metamaterials 
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Wc study the spontaneous emission of a dipole emitter imbedded into a layered metal-dielectric meta- 
material. We demonstrate ultra-high values of the Purcell factor in such structures due to a high density 
of states with hyperbolic isofrequency surfaces. We reveal that the traditional effective-medium approach 
greatly underestimates the value of the Purcell factor due to the presence of an effective nonlocality, and we 
present an analytical model which agrees well with numerical calculations. © 2011 Optical Society of America 
OCIS codes: 160.1190, 160.3918. 



Spontaneous emission rate of a ligiit source can be 
tuned by engineering its environment. A ratio of the de- 
cay rate in the media as compared to that in vacuum is 
usually termed as the Purcell factor [1], and the possibil- 
ity of a controllable change of of the radiative lifetime has 
been demonstrated in various optical systems, including 
photonic crystals and plasmonic nanostructurcs [2-4]. 

In metamaterials the Purcell factor can be greatly en- 
hanced, with two possible mechanisms of achieving high 
values of the Purcell factors. In the first case, we should 
place a dipole emitter into a metamaterial that is de- 
scribed as an uniform hyperbolic medium. Such hyper- 
bolic media are uniaxial media characterized by the per- 
mittivity tensor of the diagonal form with the princi- 
pal components being of the opposite signs [5,6]. The 
density of photonic states in such a system diverges. As 
a result, a huge Purcell factor can be reached, and its 
values being determined cither by losses and inhomo- 
geneity of the medium [6] or by a spatial extent of the 
source [7] . The second mechanism of achieving high Pur- 
cell factors employs the excitation of surface plasmon po- 
laritons (SPPs) at the metal-dielectric interfaces inside 
a metamaterial [4, 8] . 

In this Letter, we study the spontaneous emission of a 
dipole emitter imbedded into a layered metal-dielectric 
metamaterial (see Fig. 1) and reveal the possibility for 
ultrahigh values of the Purcell factor due to a high den- 
sity of states with hyperbolic isofrequency surfaces. In 
particular, we demonstrate the dramatic dependence of 
the Purcell factor near the SPP resonance on the ratio of 
the layer thicknesses and the dipole orientation, so that 
the Purcell factor can have either maximum or minimum 
at the plasmon frequency for high and low metal filling 
factors, respectively. This revealed mechanism of the en- 
hancement of the Purcell factor has an interface nature, 
and thus it can not be described within any traditional 
homogcnization approach. However, here we present an 



analytical model which agrees well with our numerical 
results. 

The hyperbolic medium at certain frequency ranges 
can be realized by layered metal-dielectric nanostruc- 
turcs. Through conventional effective-medium approach 
a metal-dielectric nanostructurc formed by pairs of lay- 
ers with permittivities ei,£2, and thicknesses di,d2 
(sec Fig. 1) can be modeled as an uniaxial anisotropic 
medium with permittivity tensor of the form: 



Eoff 




eidi + £2^2 
di + d2 
e^^di + e. 



di + d2 



At the same time, one of the key properties of such a 
medium is its strong nonlocality due to excitation of SPP 
modes at metal-dielectric interfaces [9]. 
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Fig. 1. Geometry of the layered metamaterial formed by 
alternating metal and dielectric layers. Arrows show the 
orientations of a dipole emitter (a red filled circle). 

Such structures were a subject of extensive theoreti- 
cal and experimental studies being suggested for many 
applications such as superlenses with subwavelength res- 
olution [10], a simple realization of the so-called hyper- 
lens [11] and optical nanocircuits [12], as well as invisibil- 
ity cloaks [13] . For a certain ratio of the layer thicknesses, 
one can reach the hyperbolic regime when high Purcell 
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factors can be expected [8]. In particular, it was demon- 
strated that the effective medium approach is not always 
applicable for the calculation of the radiation emission, 
and it can overestimate the Purcell factor [8]. In the 
regime we consider here the effective-medium approach 
underestimates the Purcell factor, so that the high Fur- 
cell factors can be attained even in the structures which 
arc not hyperbolic when being treated as effective media. 

We study a system shown schematically in Fig. I. A 
point dipole emitter is placed in the center of a dielectric 
layer of the periodic metal-dielectric layered nanostruc- 
ture, the simplest version of a metamaterial. In such a 
planar system, we consider independently two orienta- 
tions of the dipole: parallel (px) and perpendicular (pz) 
to the interface planes. Although the system is periodic, 
it can be formally treated as a cavity, where the central 
dielectric layer (labeled as the layer 1) with the dipole 
is surrounded by two identical mirrors. One can employ 
then the standard results of the cavity theory [14] to cal- 



culate the Furcell factors 
where 



1 + Rx, andF„^ = l + R^ 



i?^ = - Re 
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Rz = -3 Re 



where k ~ ^/eiw/c and fci^'^^ = {£i^2W^/c^ — fc^}^^^ 
Integration in Eqs. (1) is performed over the in-plane 
component of the wavevector . The amplitude reflec- 
tion coefficients of the semi-infinite periodic mirrors for 
s and p polarizations, and Vp are readily calculated 
(see, e.g. [15]) 



(la) 
(lb) 



(2) 



where ri, ti are the reflection and transmission coef- 
ficients for one period and a given polarization, D ~ 
di + d2 is the structure period, and k±_ is the Bloch 
wavevector defined by the dispersion relation: 



cos(fc_L£') 



3(4^)di)cos(A:(2)rf2)- (3) 
sin(fc(i)di)sin(fc(2)d2) 



where Zf = k'i\ Zf = Si/kz'' ,i = 1,2. In our numer- 
ical calculation we chose silver as metal with the dielec- 
tric constant described by the Drudc model [16]: 62 = 
£00 — i^p/w(a; -|- ij) with £00 ~ 4.96, ujp sa 8.98 cV,7 w 
0.018 eV, while the dielectric constant £1 was set to unity. 
The structure period is fixed to Z) = 30 nm, and only 
a ratio of the thickness of silver and dielectric layers is 
varied, rj = d2/di. 

The dependence of Furcell factors (1) on the frequency 
for different metal thicknesses is shown in Fig. 2. For bet- 
ter presentation, the Furcell factor is scaled by cube of 



Z^ 

Zl'P 



the light frequency, since the vacuum decay rate r^J is 
proportional to w'^ [17]. We observe that for the case 
when the silver layer is thicker than the dielectric layer, 
there is a sharp maximum of the Furcell factor at the 
frequency of the SFF resonance where the condition 
Re £2 = — 1 is fulfilled. However, for the case of a wider 
dielectric layer, the picture changes dramatically, and 
we observe a local minimum at the frequency of the SFF 
resonance. Next, we examine the latter case in detail. 
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Fig. 2. Furcell factor for different layer thicknesses vs. 
frequency. (a,b) Two orientations of the dipole emitter: 
parallel and perpendicular to the layers, respectively. 



The solid line in Fig. 3(a) depicts the Furcell factor 
calculated for the in-plane dipole for 77 — 0.5. These 
results reveal two characteristic effects. First, the Fur- 
cell factor can be enhanced even in the non-hyperbolic 
case. Second, the Furcell factor exhibits a sharp dip 
near the SFF resonance frequency. To explain this puz- 
zling behavior, in Fig. 3(b) we plot the isofrequency 
contours fcj_(fc||) for the close-to-resonance frequencies. 
We observe two distinct contours: an elliptic contour 
k\/ei_ + fc^/£|| = {oj/cY predicted by the effective- 
medium theory, and an additional hyperbolic contour. 
The importance of this contour has been mentioned 
in Rcf. [9], where the isofrequency contours of metal- 
dielectric nanostructures were studied by the exact and 
effective-medium approach, and they attributed to the 
effect of nonlocal response of such layered structures. 
It is this extra contour that determines a dramatic en- 
hancement of the Furcell factor in Fig. 2 in the regime 
when the structure is not hyperbolic. 

Next, we derive a closed form expression for the cor- 
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responding contribution to the Purcell factor, 
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■ in^, 5k\\ = 
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(4) 



(5) 



which is valid in the vicinity of the plasmonic resonance, 
so that ^ ^ 1. Here is the center and Sk\\ is the 
width of the hyperbohc contour. Naturahy, the Purcell 
factor is determined by the density of states (fc||''^)^5/c|| , 
and the extra factors in Eq. (4) depend on the the ori- 
entation and position of a dipole in the dielectric layer. 
The approximate expression (4), shown in Fig. 3(a) re- 
produces well the Purcell factor calculated numerically. 

The position of the center of the hyperbolic contour, 
found numerically from the equation cos[k±{k'^l^^)D] = 0, 
is plotted in Fig. 3(c), and it is well described by the 
analytical expression (4). The width of this contour as 
a function of the frequency is also reproduced by the 
function 6k\\ in Eq. (4) [see Fig. 3(d)]. Thus, this results 
uncover the origin of the dip in the Purcell factor at the 
surface plasmon frequency: while the position of the hy- 
perbolic contour increases logarithmically when ^ — oo 
as we approach the plasmon frequency, the width of the 
hyperbolic contour demonstrates a power-low decrease. 
Therefore, as a product of these two values the Purcell 
factor should vanish at the SPP resonance, but it re- 
mains finite due to losses in silver. It is also worth men- 
tioning that the effect of the enhancement of the Pur- 
cell factor demonstrated here is purely nonlocal, since 
the effective-medium approach can predict only the first 
elliptic contour [dashed line in Fig. 3(b)] and thus the 
resulting integral (1) will be small. 

In conclusion, we have studied the spontaneous emis- 
sion process and radiation rate of a dipole imbedded 
into a layered metal-dielectric metamaterial. We have 
demonstrated that the effective-medium approximation 
underestimates the radiative decay rate in such struc- 
tures which is greatly increased in comparison with a 
bulk dielectric or metal. We have predicted, both ana- 
lytically and numerically, ultrahigh values of the Purcell 
factor in such metal-dielectric structures due to a high 
density of states with hyperbolic isofrequency surfaces. 
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Fig. 3. (a) Purcell factor of a dipole oriented perpendic- 
ular to the layers for rj = 0.5. Red line corresponds to 
expression (4) ; (b) Isofrequency contours for the eigen- 
modes of the structure, Uw = 3.57 eV. (c) Position of the 
hyperbolic contour vs. frequency. Red line corresponds 
to expression (5) ; (d) Width of the hyperbohc contour 
vs. frequency. Red line corresponds to expression (5). 
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